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Abstract. Molecular gas in external galaxies is a subject of crucial im- 
portance for observational and theoretical studies of galaxy formation. 
Compact molecular gas around an active galactic nuclei (AGN) is ex- 
pected to be an energy budget of AGN and/or the nuclear starburst. 
Recent observational studies suggest that line ratios in millimeter and 
submillimeter band may be a good tool to reveal the long-standing 
question on the origin of activity - AGN or nuclear starburst. We 
have constructed a powerful "telescope" of theory, three-dimensional 
nonLTE line transfer code, preceding the high resolution and sensitiv- 
ity observations such as ALMA. 



1 Introduction 

Current progress in millimeter and submillimeter observations have found pecu- 
liar signature in ratios of molecular lines towards numbers of compact gas at the 
centers of active galaxies. The pioneering works of Kohno et al. (|20011 12005P 
showed that some of nearby Seyfert galaxies have enhanced ratio of -Rhcn/hco+ 
as high as 2.5 or more. Following observations of various molecular lines have 
given rise of discussion of peculiar chemistry such as photon dominated region 
(PDR) or X-ray dominated region (XDR) (Meijerink et al. [20051 [20071 Aalto 
et al. 120071 1 However, since the compact molecular core is not resolvable with 
current mm/submm observational instruments, one has to take into account of 
internal structures within a core. Hydrodynamic simulations present highly in- 
homogeneous density and temperature structure within a compact molecular gas 
at the center of a galaxy (Wada & Tomisaka I2005P . In an inhomogeneous molec- 
ular gas naturally expected is the inhomogeneous population distribution. We 
performed three-dimensional non-thermodynamic equilibrium (nonLTE) molecu- 
lar line transfer simulations, and examined the effects of excitation status inside a 
clumpy torus. 
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2 Ratio of HCN and HCO+ Lines in Millimeter Band 



2.1 Simulations 

We first performed hydrodynamic simulations of molecular gas in a steady grav- 
itational potential of supermassive black hole of A/smbh — lO^M© and galactic 
halo (see for details, Wada & Tomisaka !MI5l Yamada et at, I2007bp . The hy- 
drodynamic simulation explicitly includes radiative cooling, feedback from super- 
novae, strong UV heating (Go = 10) for thermal evolution as well. Resultant 
molecular ISM has a highly inhomogeneous structure in density, temperature, 
and turbulent velocity. Next we calculated three-dimensional nonLTE line trans- 
fer of HCN and HCO"*" pure rotational lines using a snapshot data of hydrody- 
namic simulation. We assumed spatially uniform molecular abundance distribu- 
tion in order to examine the effect of clumpy excitation conditions on the line ratio 
^HCN/HCO+ = ^HCN(i-o)/-fHCO+(i-o) Separately from chemistry models. 



2.2 Intensity Distribution and -Rhcn/hco+ 




Fig. 1. Integrated intensity distribution of HCN(l-O) line of line transfer simulation with 
y = 2 X 10^' (left panel). Right panel shows integrated intensity averaged over the field- 
of-view as a function of y. Line with squares denotes HCO^(l-O) line, and line with 
crosses denots HCN(l-O) line, respectively. 



In Figure [T] (a) we present integrated intensity distribution of IICN(l-O) (/hcn) 
line. The distribution of /hcn shows clumpy structure, reflecting the highly in- 
homogeneous internal structure of the model torus. The typical size of the bright 
spots in Fig. [T] is « 0.03" if the torus is located at the distance of D = 20 
Mpc. Thus though current mm-telescopes smear out these internal structures, 
high angular resolution of forthcoming telescope such as Atacama Large Millime- 
ter/submillimeter Array (ALMA) will reveal inhomogeneous structures in molecu- 
lar tori in distant galaxies. The FiglT](b) plots integrated intensities of HCN(l-O) 
and HCO"'"(1-0) avraged over the "field of view" of our simulations as a func- 
tion of fractional molecular abundance y. For the same value of y, (/hco+) 
always larger than (/hcn)j which means that a ratio of -Rhcn/hco+ ^ observed 



M. Yamada : Line Transfer Study of Extragalactic ISM 



3 



in some of nearby galaxies requires onverabundant HCN compared with HCO"*" 
(j/HCN ~ 10 X yHCO+)- This result and recent chemical studies of PDR and XDR 
do not seem to support the original idea of mm/submm molecualr line diagnos- 
tics of AGN/starbust connection (e.g. Kohno et al. I200ip . Other models, such 
as mid-InfraRed (MIR) photon pumping (e.g., Gracia-Carpio et a^ l2006l Aalto 
et a^ l2007p or hot core chemistry (e.g., Stauber et a/ 12005p have been proposed, 
though, theoretical model description of molecular line ratios been controversial. 



3 Multi-Transition Analysis in Inhomogeneous Torus 

Current progress in submm observations enables multi-transion analysis and probe 
of warm and dense gas at the centers of external galaxies. In Figure [2] ratio of 
integrated intensities (^hcn(4-3)) and (^hcn(i-o)) (hereafter denoted as -R43/10) 
is plotted as a function of fractional molecular abundance y. In our simulation 
results ratio -R43/10 clearly decreases with y. 
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Fig. 2. -R43/10 as a function of viewing angle with respect to rotational axis of the torus. 
Right panel shows population ratio of J = 4 and 1, or excitation temperature of rcx(41). 



On the other hand, one-zone analysis predicts a reverse trend : if ?; is small. 



optical thickness tq is expected to be small, then the ratio R 



■43/10 



will be 



R. 



43/10 
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0.3, 



(3.1) 



where 7oj is coUisional excitation constant from energy level J = to J. If y is 
large, tq will be also large, and i?43/io approaches unity. 

This result indicates the importance of the inhomogeneity in structure of the 
torus. We model the clumpy torus as two-phase ISM consists of 1) dense clumps 
with n ~ ricrit and 2) tenuous ambient medium with n <^ «crit- If we assume 
high kinetic temperature (/cB^kin ^ hv) and optically thin over a whole region. 
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intensity ratios become (see eqs. [9] and [10] of Yamada et al. I2007ap 




(3.2) 




(^1^43 E,7>4 70,7 ?lO> 



(3.3) 




Mean ratio (i?43/io)an in a- synthetic observation beam can be written as a kind 
of average of (i?43/io>d and (i?43/io>t: 



where rj denotes a volume filling factor of dense clumps. Equation (|3.4p means 
that if rj is sufficiently smaller than unity, i?43/io can take a value of order unity 
even if the medium is optically thin. In our radiative transfer simulations, "dense 
clumps" can safely be defined as regions with n > 10^ cm~'^ (close to Ticrit of 
HCN(l-O) transition), and ?7sim is indeed « 10~^. Furthermore, nonLTE pop- 
ulation (occasional tq < for J =1-0 transitions) and inhomogeneous optical 
thickness distribution comphcate the -R43/10 distributions (see right panel of Fig. 
[2|). Our analysis demonstrates the importance of inhomogeneity of molecular gas 
in estimation of physical properties from line ratios from unresolved observational 
data. Inherent inhomogeneity should be paid much attention when interpreting 
line ratio in terms of chemical evolution of molecular gas. 

4 Summary 

We examine molecular line diagnostics of molecular gas in the center of exter- 
nal galaxies using three dimensional radiative transfer simulations. Our results 
demonstrate that realistic inhomogeneous molecular gas can produce line ratio 
values unexpected from one-zone analysis, and emphasize the ability and appli- 
cability of our approach, that is, compilation of high-resolution hydrodynamic 
simulations and three-dimensional radiative transfer calculations ("telescope of 
theory") towards forthcoming ALMA era. 
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